Quantum fluctuations and condensate fraction during the time-of-flight expansion 
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The quantum fluctuation effects in the time-of-flight (TOF) experiment for a condensate released 
from an optical lattice potential is studied within the truncated Wigner approximation. By investi- 
gating both the spatial and momentum density distributions, we find that the condensate fraction 
decreases monotonically in time and hence cannot be measured in the standard TOF image. We 
then propose a semi-quantitative analysis for such dynamical quantum depletion process. Our study 
shows a universal algebraic decay of the true condensate fraction, and have a very good agreement 
with numerical results. We also discuss possible methods to determine the condensate fraction inside 
the optical lattice, and its implication to the TOF experiments in higher dimensional systems. 
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Experimental realization of the superfluid to Mott in- 
sulator transition in ultracold bosonic atoms opens a new 
horizon of strongly correlated physics [l] and Q . Among 
various kinds of experimental techniques, time-of-flight 
(TOF) experiment is the most important tool to observe 
possible quantum orders and related correlations 
It is generally believed that the TOF image can be inter- 
preted as the momentum distribution of the atom cloud, 
if (i) the final size of the atom cloud is much larger than 
the initial one and (ii) the interaction effect during the 
expansion is negligible. As a result, the sharp interfer- 
ence peaks near the zero and Bragg wavevectors can be 
interpreted as a signature of a Bose-Einstein condensa- 
tion inside the optical lattice. 

Recently, Gerbier et ai. |6J show that the condition (i) 
is usually not satisfied, because the required time scale 
is much larger than the feasible value in the present ex- 
periments. As for the condition (ii), it is argued that the 
interaction effects are negligible (or at most in the mean- 
field level) due to the sudden drop of the particle density 
Different from the condensate in a single trapping 
potential, however, particles loaded in an optical lattice 
contain a much stronger interaction energy (and hence 
quantum fluctuation). In principle, the interaction effect 
can be removed by using Feshbach resonance, but such a 
far-from-equilibrium quantum dynamical problem itself 
is still very interesting and worthy of further investiga- 
tion. It may also provide valuable insights for studying 
other quantum orders through the TOF experiment. 

In this paper, we numerically and analytically in- 
vestigate the quantum dynamics of the TOF experi- 
ment for an interacting condensate released from a one- 
dimensional (ID) optical lattice potential. The quantum 
fluctuations are included within the truncated Wigner 
approximation (TWA). The common adapted empirical 
definition of the condensate fraction, measured from the 
bi-modal fitting of a spatial distribution, is found to sat- 
urate at a value substantially smaller than the free ex- 
pansion result. On the other hand, the true condensate 
fraction, which can be unambiguously determined in the 



momentum space at any moment, decreases monotoni- 
cally in time due to the dynamical quantum depletion. 
We find that the decaying condensate fraction can be fit- 
ted very well by a simple algebraic function with a uni- 
versal exponent, showing an interaction induced quan- 
tum decoherence effects. We also show the 2D results 
and find a similar quantum depletion phenomenon. Our 
results provide a deeper understanding about such a far- 
from-equilibrium quantum dynamical process and can be 
applied to studying other quantum orders. 

We consider 87 Rb atoms loaded in a ID optical lat- 
tice with the scattering length, a s — 5.45 nm, the lattice 
constant, a — 0.5 fim, and the lattice strength, Vb = 10 
En, where Er = 2irh x 2.3 kHz is the recoil energy. The 
transverse confinement radius is R = 0.11 /im, the on- 
site interaction energy is U = 0.2£^. When t > 0, the 
lattice potential is turned off and atoms start to expand 
with the same transverse confinement (with an effective 

ID interaction strength, g\ = ^FiPm ' wnere m is the 
atomic mass). We use three approximations for numer- 
ical studies: non-interacting expansion, Gross-Pitaevskii 
equation (GPE), and the truncated Wigner approxima- 
tion (TWA), respectively. Note that the TWA can be 
shown to be exact in the short-time limit as well as in 
the weak interaction limit [9l-flTj|. and therefore is very 
appropriate for the study of TOF dynamics. In all of 
our following calculations, we apply the same initial con- 
dition, a pure condensate state with ten particles in ten 
lattice site, for simplicity. Application of our results to 
different initial states (say Mott state or with finite ther- 
mal/quantum depletion) is feasible, but is not our in- 
terest in this paper. Finally, we discretize the optical 
lattice potential by eight mini-cells in space, and calcu- 
late the expansion dynamics of the whole system up to 
2048 mini-cells. In other words, for the initial condensate 
of ten lattice sites, the atom cloud is allowed to expand 
at most 25 times larger than its initial size within 100 
h/Ejj ~ 7 ms. 

In Fig. QIal)-(a4), we show the calculated spatial den- 
sity distributions at different expansion time. We note 
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FIG. 1. (Color online) Time evolution of the spatial (al)- 
(a4) and momentum (bl)-(b4) distributions for an interacting 
condensate released from a ID optical lattice. The expansion 
time (t) is in the unit of h/Eu. The dashed black, dotted red, 
and solid blue lines are the non-interacting, GPE and TWA 
results, respectively. The dash-dotted green lines in (a3) and 
(a4) are proportional to the square of the expanding Wannier 
function. Other parameters are described in the text. 



that both the GPE (dotted red lines) and the TWA 
(solid blue lines) results show a strongly broadened cen- 
tral peak compared to the non-interacting result (dashed 
black lines), while the difference becomes much weaker 
in the two Bragg peaks. To quantify such an interaction 
effect, we define the following renormalized peak ratio, rj: 
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where w(k) is the Fourier transform of the Wannier func- 
tion at the wavevector k,x(k) = htk/m , and p{x,t) 
is the particle density at position x at time t. We 
note that r] is exactly equal to one for a free expan- 
sion. This can be shown from the free expansion result: 
p(x(h),t) oc |w(k)| 2 iS(k) |6j, where the structure func- 
tion, <S(k), is a periodic function with the period 2n/a, 
i.e., 5(0) = S(2it/a). Nevertheless, when the interaction 
is considered, our TWA result shows that n can be re- 
duced to about 0.65 at t = 75H/ER for the simulations 
in Fig. Q] In fact, similar experimental result has been 



observed in Ref. [12[ , where the population ratio of the 
first Bragg peak to the central peak is 35% larger than 
the ratio of a free expansion (i.e., r/(t) — ~ 0.74) at 
t = 29.5 ms. A systematic study on r/(t) should provide 
more information about the meanfield interaction effects 
during the TOF experiment. 

At the first glance, the similar real space distributions 
of the GPE and TWA in Fig. QIal)-(a4) seem to sug- 
gest that the quantum fluctuation is negligible. However, 
from the momentum distribution shown in Fig. [TJbl)- 
(b4), one can see that the TWA result deviates from the 
non-interacting or the GPE results seriously: the con- 
densate particles in the central peak are depleted and 
more non-condensate particles are generated with finite 
momenta, filling up the valleys of the initial momentum 
distribution. Note that the GPE result also has valleys 
after a longer time expansion, which simply results from 
the finite size effect of the initial wavcfunction. In order 
to quantitatively investigate such a dynamical quantum 
depletion process, here we define the condensate fraction 
f c {t) to be the ratio between the number of condensate 
particles N c (t) near the zero momentum at time t to its 
initial value, i.e., f c (t) = N c (t)/N c (0). We use two in- 
dependent numerical methods to extract f c (t) from our 
TWA results: one is from the standard bi-modal fitting 
of the central peak in the momentum distribution (see 
Ref. [IH and the inset of Fig. [5]), and the other one is by 
calculating the probability to find particles in the con- 
densate coherent state [14]], which is supposed to follow 
the GPE result if no quantum fluctuation is included. 
In Fig. HJa) we show the calculated results obtained 
from these two methods in the dotted red and solid blue 
lines, respectively, and find a very good agreement: the 
condensate fraction decreases monotonically during the 
TOF expansion and the decay rate becomes larger for 
a stronger interaction. This quantum fluctuation effect 
cannot be obtained from our GPE (mean-field) results at 
all. We also note that when this method will fail auto- 
matically if the initial state is a normal or Mott state so 
that iV c (G) = 0. 

For comparison, now we further calculate the pseudo- 
condensate-fraction, f c (t), which is extracted from the 
bi-modal fitting in the spatial density distribution (see 
Ref. [l3|), as the experimental convention. In Fig. [2ja), 
we show the calculated f c (t) (raising green lines) for three 
different interaction strengths: it grows monotonically in 
time and eventually saturates at a value, / c (oo), which 
is equal to one for a free expansion but smaller than one 
with a finite interaction. This result indicates that if the 
expansion is non-interacting, one can correctly measure 
the true condensate fraction after a long time expansion 
(/c(oo) = f). On the other hand, when the interaction is 
considered, f c {t) saturates at a value / c (oo) < f , which 
cannot represent the true condensate fraction inside the 
optical lattice. Since the results shown in Fig. [Ha) are all 
calculated from the same initial condition, we conclude 
that the spatial density distribution (i.e., the TOF im- 
age) cannot give a good measurement of the condensate 
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FIG. 2. (Color online) (a) Decay of the condensate fraction, 
f c (t) for different interaction energies, U/Er = 0.1 and 0.3. 
The dotted red and solid blue lines are obtained from the 
bi-modal fitting p]§| in momentum distribution and wave- 
function overlapping [Jf] , respectively, while the dashed black 
lines is fitted by Eq. (|3]) . The three raising green lines are the 
pseudo-condensate-fraction, f c {t), for U/Er = 0, 0.1, and 0.3 
(in the order from top to bottom). Inset is the momentum 
distribution with the bi-modal fitting cures, (b) Same as the 
case in (a) for U /Er = 0.3, but the interaction is turned off 
at t, = hh/ER (upper line) and t* = 40H/Er (lower line), 
respectively. 



fraction during the TOF. 

Although the time-dependence of the true condensate 
fraction, f c (t), and the pseudo-condensate fraction, f c (t), 
are quiet different, their close relationship can be revealed 
by switching off the particle interaction immediately at 
time t* during the expansion. As shown in Fig. EJb), 
f c (t) evolves eventually to a value set by the true con- 
densate fraction at f*, i.e., / c (oo) = fc{t*), no matter its 
value is larger or smaller than f c {t*)- This is simply be- 
cause, after i», particles expand freely in space and there- 
fore the momentum distribution at is transformed to 
the spatial distribution. This also confirms the decay of 
the condensate fraction. Note that, without turning off 
the interaction during the expansion, the true condensate 
fraction, f c (t), is hidden in the momentum distribution 
and therefore cannot be directly observed in the TOF 
image. However, turning off the inter-particle interac- 
tion can also reduce the expansion speed, making it more 
difficult to reach the desired saturated value, / c (oo). 

The quantum depletion during the TOF can be viewed 
as a radiative process in the momentum space, where 
a pair of condensate particles are scattered out of the 



condensate and transfer their interaction energy toward 
the kinetic energy. Such dynamical process is gov- 
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where 



go is the interaction strength in the D-dimensional 
free space, 6 q is the boson operator, and Q d is 
the system volume. Here the scattering of parti- 
cles in the Bragg peaks are neglected for simplicity. 
Denoting the many-body wavefunction at time t by 
|\&o(*))i the many-body wavefunction at time t + dt 
(dt — > + ) can be described by the following form: 
|*o(< + dt)) = B (dt)|*o(t)> + E q Bq(dt)l*q(*))i where 
|* q (i)> = 7V c (i)-i& q ^ q 6 6 |>F (i)) and B^dt) is the as- 
sociate coefficient. The non-condensate particle num- 
ber, (\I r o|& q & q |\E r o), has been neglected when compared 

to N c (t) = (<Fo(t)|&Jfe | , I'o(£)) in the normalization con- 
stant. Therefore, the condensate particle number at t+dt 
becomes N c (t + dt) = N c (t) - 2 £ q |-B q (cft)| 2 . Giving 
Bi(0) = 5i t o, we can apply the spirit of Fermi-Golden 
rule for dt — > 0, and obtain the decay rate of the conden- 
sate density n c (t) = N c (t)/flD by 
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where Cjj oc A (jp) D ^ 2 is a constant, and fijj is the chem- 
ical potential in the D-dimensional space. We note that 
in general /i d depends on the interaction strength, global 
trapping potential, and dimensionality etc. For example, 
in a uniform 3D system, fi ~ gzn c in the weakly interact- 
ing limit, while there is no condensate particles in a uni- 
form ID system. For a ID finite size interacting system 
here, without loss of generality, we may simply assume 
flu = erjn c , where to should depend on the interaction 
strength as well as system size. After integrating over 
this differential equation, we obtain the condensate frac- 
tion, f c (t) = n c (t)/n c (0), as following: 
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with a non- universal fitting parameter here Co = 
Cd^q^ 2 1 (depending on finite size effect, finite temper- 
ature, initial non-condensate particles etc.). It is impor- 
tant to note that, instead of an exponential decay ex- 
pected in a regular radiative process, the many-body in- 
teraction matrix element, H' I: couples the decay rate to 
the condensate density, so that it becomes an algebraic 
decay instead. If the initial state is not a pure conden- 
sate, Eq. ([3]) should be still valid if only away from the 
regime of quantum phase transition or Mott phase. Us- 
ing Eq. © to fit our TWA results shown in Fig. rjfa) 
with a fitting parameter Cr>, we find a very good agree- 
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FIG. 3. (Color online) Momentum distribution of a TOF 
expansion in 2D space along y — at t = H/Er. The dashed 
black, dotted red, and solid blue lines are the non-interacting, 
GPE and TWA results, respectively. (The first two cannot 
be distinguished here.) The two insets show the momentum 
distributions for TWA at t = and t = H/Er, respectively. 
Here the interaction energy is larger (U — 0.45-Ek) in order 
to demonstrate the quantum depletion effects within a short 
time expansion. 



merit to confirm the quantum fluctuation effects on the 
condensate fraction during the TOF experiment. 

As for higher dimensions, in Fig. [3J we also show 
the momentum distribution of atom clouds expanding 
in a 2D system from an initial 11 x If lattice sites (to- 
tal 128 x 128 mini-cellls) . Due to the limitation of our 
computational resources, only results of a very short time 
expansion can be obtained before atoms reach the bound- 
ary. No quantitative information can be extracted yet. 
However, by comparing the momentum distributions in 



the GPE and TWA results, it is easy to find the decay of 
condensate density even in such a short-time regime due 
to the quantum fluctuations. 

Before concluding, we mark several aspects to discuss: 
(i) We have also calculated the TOF dynamics of twenty 
atoms in ten lattice sites (not shown here), and find a 
stronger decay of f c (t) for the same interaction strength. 
However, results of twenty particles in twenty wells are 
found almost the same as results of ten particles in ten 
wells. This shows that the system size is not relevant 
here, and its dependence of the particle density (rather 
than total particle number) is also consistent with the 
nature of quantum fluctuations, (ii) Since TOF is known 
as a far-from-equilibrium dynamical process, the integra- 
bility of ID Bose gas should be irrelevant to any of our 
results here. This is also confirmed by the similar de- 
pletion process in 2D systems shown in Fig. |3j (iii) 
From Eq. ((3J, it is easy to see that for an expansion in 
higher dimensional systems, the exponent of decay be- 
comes smaller, due to the smaller density of states in the 
long- wavelength limit. This is therefore also consistent 
with the general principle that the quantum fluctuation 
effects become smaller in higher dimensional systems. 

In conclusion, we numerically and analytically study 
the interaction effects for an expanding condensate re- 
leased from an optical lattice. The condensate fraction 
has an algebraic decay in time due to the quantum fluctu- 
ation effects. Our results provide a deeper understanding 
to the quantum fluctuation effects during the time-of- 
flight experiment. 

We appreciate the discussion with T.-L. Ho, C. Chin, 
I. Spielman, T. Proto, and G. Modugno. We also thank 
P.-C. Chen for his kind assistance on the computation. 
This work is supported by NSC(Taiwan). 



[I] M. Greiner, et al., Nature 415, 39 (2002); F. Gerbier, 
et al, Phys. Rev. Lett. 95, 050404 (2005); I. B. Spiel- 
man, W. D. Phillips, and J. V. Porto, Phys. Rev. Lett. 
98, 080404 (2007); V. Guarrera, et al., Phys. Rev. Lett. 
100, 250403 (2008); G. Roati, et al, Phys. Rev. Lett. 99, 
010403 (2007). 

[2] I. Bloch, J. Dalibard, and W. Zwerger, Rev. Mod. Phys. 

80, 885 (2008). 
[3] I. B. Spielman, W. D. Phillips, and J. V. Porto, Phys. 

Rev. Lett. 100, 120402(2008). 
[4] R. B. Diener, et al, Phys. Rev. Lett., 98, 180404 (2007); 

T.-L. Ho and Q. Zhou, Phys. Rev. Lett. 99, 120404 

(2007). 

[5] M. W. Zwierlein, et al, Nature 435, 1047(2005); V. 

Schweikhard, et al., Phys. Rev. Lett. 93, 210403(2004); 

J.K. Chin, et al, Nature 443, 961 (2006). 
[6] F. Gerbier, et al., Phys. Rev. Lett. 101, 155303(2008). 
[7] G.-D. Lin, W. Zhang, and L.-M. Duan, Phys. Rev. A 77, 

043626(2008); W. Yi, G.-D. Lin, and L.-M. Duan, Phys. 

Rev. A 76, 031602(R) (2007). J. N. Kupferschmidt, and 

E. J. Mueller, I arXiv: 1006 .20421 (unpublished) 
[8] Y. Castin and R. Dum, Phys. Rev. Lett. 77, 5315 (1996); 



Yu. Kagan, E.L. Surkov, and G.V. Shlyapnikov, Phys. 
Rev. A 55, R18 (1997). 
[9] D.F. Walls and G.J. Milburn, Quantum Optics (Springer- 
Verlag Berlin Heidelberg New York, 1994). 
[10] M.J. Steel, et al, Phys. Rev. A 58, 4824 (1998); A. Sina- 
tra, et al., J. Mod. Opt. 47, 2629 (2000); Phys. Rev. Lett. 
87,210404 (2001); J. Phys. B: At. Mol. Opt. Phys. 35, 
3599 (2002). 

[11] A. Polkovnikov, Phys. Rev. A, 68, 033609(2003); Phys. 
Rev. A, 68, 053604 (2003); A. Polkovnikov and D.-W. 
Wang, Phys. Rev. Lett. 93, 070401 (2004). 

[12] P. Pedri, et al, Phys. Rev. Lett. 87, 220401 (2001). 

[13] Our bi-modal fitting process is following: We first set a 
reasonable momentum fco and fit the distribution by a 
Gaussian function for |fc| > fco. After substracting this 
first fitting curve, we fit the residual distribution by a 
second Gaussian function for jfc| < fco. The area of the 
second fitting result gives N c (t), which is checked to be 
independent of fco. Condensate particles at Bragg peaks 
are neglected. Similar method is applied to extract f c {t) 
from the real space distribution. 

[14] The condensate density obtained from the wavefuction 



■5 



overlapping method is to calculate (6j p (£)6 G p(t)), where ator. See also M.J. Steel, et al, Phys. Rev. A. 58, 4824 

bGp(t) = J ipGp( x !~t)iP( x )dx is the operator projected (1998). 
onto a GP state, ipGp(x,t). Here i>(x) is the field oper- 



